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There is a need for covariant solutions of bound state equations in order to construct realistic QCD 
based models of mesons and baryons. Furthermore, we ideally need to know the structure of these 
bound states in all kinematical regimes, which makes a direct solution in Minkowski space (without 
any 3-dimensional reductions) desirable. The Bethe-Salpeter equation (BSE) for bound states in 
scalar theories is reformulated and solved for arbitrary scattering kernels in terms of a generalized 
spectral representation directly in Minkowski space (K. Kusaka et al., PRD 51, 7026 '95). This 
differs from the conventional Euclidean approach, where the BSE can only be solved in ladder 
approximation after a Wick rotation. 
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1 Introduction 



The Bethe-Salpeter Equation (BSE)B describes the 2-body component of bound-state struc- 
ture relativistically and in the language of Quantum Field Theory. It has applications in, for 
f-» example, calculation of electromagnetic form factors of 2-body bound states and in studies 

of relativistic 2-body bound state spectra and wavefunctions. 

BSEs have been solved previously for separable kernels and for ladder scattering kernels. 
Solutions for BSEs have also been obtained for QCD-based models of meson structure in 
Euclidean space; these solutions must be analytically continued to Minkowski space. If we 
4h ' solve the BSE with a non-ladder scattering kernel or dressed propagators for the constituent 

^p" 1 ' bodies, the validity of this procedure (known as the Wick rotation!!) is highly non-trivial, 

Qh! and so direct solution in Minkowski space without unnecessary approximations is preferable. 

Here we outline such a method for scalar theories, based on the Perturbation Theoretic 
Integral Representation (PTIR) of Nakanishiu. 

The PTIR is a generalisation of the spectral representation for 2-point Green's functions 
to n-point functions; for a particular renormalised n-point function, the PTIR is an integral 
representation of the corresponding infinite sum of Feynman graphs with n external legs. 

Our approach involves using the PTIR to transform the equation for the proper bound- 
state vertex, which is an integral equation involving complex distributions, into a real 
integral equation. This equation may then be solved numerically for an arbitrary scattering 
kernelQ. We have to date applied this formalism to the case of a ladder kernel, as well as the 
so-called "dressed ladder" kernel in which we add self-energy corrections to the propagator 
of the exchanged particle. These results have been verified by direct comparison to those 
obtained in Euclidean space by previous authors El. 

An example of one such scalar theory to which our formalism may be applied is the (j) 2 a 
model, which has a Lagrangian density 



where g is the 4>-a coupling constant. 
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2 Formalism and PTIR 



The Bethe-Salpeter equation in momentum space for a scalar-scalar bound state with scalar 
exchange is 

*(p,P) = -D( Pl )D( P2 ) J A^^( q ,P)K(p,q;P), (2) 

where <3? is the Bethe-Salpeter (BS) amplitude, and where K is the scattering kernel, which 
contains information about the interactions between the constituents of the bound state. 
We may also write this in terms of the bound state vertex V as 

rfa-P) = / ■^D(p 1 )D(p 2 )K(p,q;P)T(q,P), (3) 

In Eq. (|2|), pi is the four- momentum of the i th constituent. We also define p = i]2Pi—r)ip2, 
which is the relative four-momentum of the two constituents, and P = p\ + P2 is the total 
four-momentum of the bound-state. The real positive numbers iji are arbitrary, with the 
only constraint being that rji + 7/2 = 1. We will use 771 = 1/2 = 772, which is a convenient 
choice when the constituents have equal mass. 

In order to convert the BSE into a real integral equation, we will need to use the PTIR 
for both the proper bound-state vertex and the scattering kernel!. The bound-state vertex 
may be represented as 

F(q, P) = r da I" dz .,„ • W 



J-l [a + m 2 - {q 2 + zq ■ P + - ie] 



The weight function p n (a, z) of the vertex has support only for a finite region of the space 
spanned by the parameters a and z. 

The representation for the vertex in Eq. is for s-wave bound states. It is straight- 
forward to generalise our arguments to higher partial waves QO. 

We have introduced a dummy parameter n, which will be of use in our numerical work 
since larger values of n produce smoother weight functions. The fact that n is arbitrary can 
be seen by integrating by parts with respect to a; in this way weight functions for different 
values of n may be connected i. 

We may use the PTIR for the bound-state vertex to derive the PTIR for the Bethe- 
Salpeter amplitude, since the two are related via 

fc(p, P) = iD(| + p)iT(p, P)iD(^ - p). (5) 

We proceed by absorbing the two free scalar propagators into the expression for the vertex, 
Eq. (||), using Feynman parametrisation. After some algebra we obtain for the BS amplitude 

*(p,P) = -• I' dz r da- (6) 
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m 2 + a — {p 2 + zp ■ P + — ie 
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To include the most general form of the scattering kernel in our derivation, we use the 
PTIR for the kernel: 



roc p 

-iK(p,q;P) = £ / d ^ / # 

pch{j,£) 

7 - Kh<7 2 + b ch p ■ q + c ch p 2 + d ch P 2 + e ch q ■ P + f ch p ■ P) - ie 

(7) 

where the kernel parameters {a c h> • • • > /ch} are linear combinations of the Feynman param- 
eters {£1,... >^6}) an d Q denotes the region of integration which is that imposed by the 
constraint Ya=iCi = 1- Here we have defined, similarly to before, q = {qi — ^2) /2. The 
expression for the kernel contains a sum over three different channels, labelled by st, tu and 
us. 



3 Derivation of Equations for Scalar Models 

We use the PTIR form of the kernel and vertex in Eq. (||), along with Feynman parametri- 
sation, to derive a real, two-dimensional integral equation for the vertex weight function 
p(a,z). We use bare constituent propagators, and using the PTIR uniqueness theorem!, 
obtain the form of the equation that can be solved numerically, 

lp n (a,z) f°° f 1 y,-- sPn(a,z) 

- — — — =/ da dz lC n (a,z;a,z) — (8) 

A a n Jo J-i a n 

We have defined here an eigenvalue A = g 2 / (47r) 2 , where g is the <f>-a coupling. We solve 
Eq. (H) numerically by iteration to convergence for the modified weight function p n (a, z)/a n . 
For the sake of brevity, we will not write down the explicit form of the kernel function /C 
here. Note that /C should not be confused with the scattering kernel, K, discussed earlier. 

As a test of our implementation of the formalism presented here, we have solved Eq. (|8|) 
with n = 2 for a ladder kernel and for a dressed ladder kernel where we add a one-loop 
self-energy correction to the propagator of the exchange (a) particle. Our eigenvalues 
agree with those obtained by Linden and Mitter @ in Euclidean space to typically 1 part 
in 10 or better, for moderate numbers of grid points. We plot the eigenvalue spectrum 
for ladder exchange for various values of the exchange mass \i in Fig. [l|. Fig. ^ illustrates 
the development of the modified weight function p2{a, z)/a 2 obtained by solving the ladder 
BSE (with n/m = 0.5) as we increase the value of the bound state mass squared, P 2 , 
towards the instability threshold where P 2 = Am 2 . As P 2 increases towards threshold, 
the binding becomes weaker, i.e., g decreases. Finally, in Fig. ^ we compare the weight 
functions for the ladder and dressed ladder cases (here the pole in the cr-propagator is 
located at ji/m = 1.0) when the binding is quite weak {^P 2 /Am 2 = 0.99). We found that 
the one- loop correction enhances the binding very slightly fby a few percent or smaller), 
and does not have a significant effect on the weight functional, as is borne out by Fig. ||. 
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Eigenvalue Spectrum 
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Figure 1: Eigenvalue spectrum for the scalar BSE. Here fj,/m is the ratio of the 
exchange mass to the constituent mass, and rj is a parameter which describes 
the "fraction of binding": r\ = (P 2 /4m 2 ) 2 , which in the rest frame of the 
bound state reduces to P°/2m. 



4 Conclusions and Outlook 

We have formulated the BSE in Minkowski space for a completely general scattering kernel. 
This has been solved numerically for the ladder and dressed ladder kernels; both sets of 
results agree to high accuracy with results using more traditional Euclidean space methods. 

We are currently implementing a kernel code which will enable us to numerically solve 
the scalar BSE for a completely general scalar kernel. As a first application of this code, 
we will solve for a "generalised ladder" kernel Q. The Wick rotation is invalid in this case 
and so solution of the BSE in Minkowski space is essential in this instance. 

To study problems of more general interest, it will be essential to extend the PTIR to 
theories containing fermions, and to find a way of incorporating confinement and derivative 
coupling into our approach. This will enable us to carry out studies of mesons in QCD, for 
example. 
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Figure 2: Development of the weight function as the bound state mass increases towards 
threshold. The left column contains the cases ??=0,0.2 and 0.6, while the right hand column 
contains the cases »?=0.8,0.95 and 0.99. 
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